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Abstract 



A self-consistent field theory study of lock and key type interactions between sterically stabilized 
colloids in polymer solution is performed. Both the key particle and the lock cavity are assumed 
to have cylindrical shape, and their surfaces are uniformly grafted with polymer chains. The lock- 
key potential of mean force is computed for various model parameters, such as length of free and 
grafted chains, lock and key size matching, free chain volume fraction, grafting density, and various 
enthalpic interactions present in the system. The lock-key interaction is found to be highly tunable, 
which is important in the rapidly developing field of particle self-assembly. 
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I. INTRODUCTION 



Self-assembly of colloidal particles into desired structures, which is driven by selective and directional interparticle 
interactions, can provide a promising route for fabrication of novel materials .11 [One particularly promising approach 
in this research area is based on utilizing "lock and key" colloidal systems. 2, 3j As such, these model systems have 
recently received substantial attention both experimentally @, Q and theoretically. 

In a recent experimental study, [2J an efficient method has been developed to produce colloidal lock particles 
containing a spherically shaped cavity. The interaction between these lock particles and complementary spherical key 
particles is comprised of two major contributions: Coulomb repulsion arising due to charge stabilization, and depletion 
attraction[9( due to the presence of the polymeric depletants in solution. The depletion attraction is the strongest when 
the radius of the key particle exactly matches that of the cavity due to the fact that the overlap of excluded volumes 
associated with lock and key particles is maximized in this case. As a result, this key-lock binding is highly specific, 
which is one of the crucial requirements for developing a successful self-assembly process. In addition to geometric 
considerations, the strength of depletion attraction can be controlled either by varying the depletant concentration 
or by adjusting temperature, thereby changing the polymeric depletant size.[|[ Hence, the key-lock interaction is 
not only selective, but also reversible and tunable, which is equally important for controlling self-assembly. In this 
regard, it was found that the binding-unbinding transition in the key-lock system is significantly sharper when steric 
stabilization (arising due to the presence of grafted polymer layers on both key and lock surfaces) is used instead of 
charge stabilization. Q The sharpness of the binding-unbinding transition is determined by considering the fraction 
of occupied lock cavities as a function of depletant concentration. Experimental data show that for charge stabilized 
key-lock systems this function grows gradually, while for sterically stabilized systems it approaches a step function, 
i.e. the binding- unbinding transition becomes much sharper. [2] 

On the theoretical side, key-lock interactions in colloidal systems have been studied using integral equation 
theory, Q, density functional theory, Q and molecular simulation techniques. Q Very recently, a highly efficient hybrid 
approach to study key-lock model systems (among other multidimensional problems) has been developed, which 
combines Monte Carlo simulation for obtaining microscopic configurations of the depletant particles and density 
functional theory for calculating the free energy. Q All these previous studies employed microscopic models based 
exclusively on hard-sphere excluded volume interactions, whereby the depletion attraction between lock and key 
particles is driven by the entropic solvent contribution. In agreement with experimental observations and geometric 
arguments based on excluded volume overlaps, it has been found that the strongest binding interaction occurs for the 
systems where the size of the key exactly matches the size of the lock cavity. However, it was also found that a simple 
Asakura-Oosawa type treatment Q based exclusively on considering excluded volume overlaps is not always sufficient. 
In particular, these studies have highlighted the importance of employing a realistic microscopic model for the solvent 
(depletant) by showing that the Asakura-Oosawa approach, which assumes ideality of the solvent, can significantly 
underestimate the selectivity of key-lock interactions. [J] 

While earlier theoretical studies of lock and key model systems have considered several different geometric shapes 
of key particles and lock cavities, the depletant has always been treated as a simple hard sphere solvent, thereby 
limiting the key-lock interaction to purely entropic depletion attraction. At the same time, in the experimental work 
the strength of key-lock binding has been tuned by changing the effective size of polymeric depletant via temperature 
changes, which indicates the potential importance of enthalpic effects in these systems. @ Furthermore, the fact that 
the sharpness of binding-unbinding transition can be increased in the presence of steric stabilization, Q points to the 
necessity of developing a microscopic model of lock and key systems that would explicitly include grafted chains on 
the surfaces of lock and key particles. The central goal of the present work is to develop such a model. 

In our earlier work, we have employed mean-field type theoretical methods (density functional theory and self- 
consistent field (SCF) theory[Toj) to study interactions between sterically stabilized colloidal particles in polymer 
melts and solutions. (ill [l2[ By comparing theoretical results with computer simulations, it was shown that theory 
yields accurate density profiles of both free and grafted chains, [12j as well as potentials of mean force (PMF) between 
sterically stabilized colloids. ll| In the present work, we employ mean-field techniques to study interactions between 
sterically stabilized lock and key particles in the presence of polymeric depletants. We compute the corresponding 
PMFs and study their dependence on the geometric matching between the sizes of key particle and lock cavity, which 
is the crucial parameter governing the strength of lock-key interactions. In addition, we study the dependence of 
PMF on several other model parameters, which are known to play an important role in governing the microstructure 
of sterically stabilized colloids in polymer melts and solutions, such as the ratio of free and grafted chain lengths. fll. [iH- 
[l6j grafting density, 0, [l3[ and enthalpic interactions of free chains with grafted chains [l2j and with the surface of 
lock and key particles. 

The remainder of the paper is organized as follows. In Section [TT] we introduce our microscopic model and 
describe SCF-based approach for calculating PMF between sterically stabilized lock and key particles. The results of 
our calculations of PMFs for a wide range of model parameters are presented in Section IIIII Section IIVI concludes the 
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Figure 1 
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FIG. 1. A schematic side-view of the microscopic model. The key particle (shown in red) is a cylinder of height H and radius 
Rk\ the lock cavity (shown in black) is a cylinder of height D and radius Ri. Both particles are uniformly grafted with polymer 
chains of length N (shown in blue) at grafting density a g and immersed in a solution of polymer chains of length P (shown in 
green) with bulk volume fraction (fit- 



paper. 



II. MICROSCOPIC MODEL AND THEORY 



Previous studies of lock-key interactions have considered various geometries for the particles, including spherical 
and ellipsoidal shapes. 4 8] In the present work, we assume cylindrical shape both for the key particle and for the 
lock cavity. In particular, the key particle is assumed to be a cylinder of radius and height H, while the lock 
cavity is characterized by radius Ri and depth D. The conditions Rk = Ri and H — D correspond to the perfect size 
matching between lock and key. The particles are immersed in a solution of polymer chains of length P with bulk 
volume fraction 4>b. In addition, the surfaces of the particles are uniformly grafted with chains of length N at grafting 
density o~ g . In what follows, we consider only coaxial arrangements of the cylinders describing key particle and lock 
cavity, both of which are aligned along the z axis. As a result, the lock-key PMF is a function of a single coordinate 
z. A schematic side-view of our microscopic model is presented in Fig. [1] 

We use SCF theory to obtain the PMF between sterically stabilized lock and key particles in a polymer solution. 
The central quantity in the SCF approach is mean-field free energy, which is expressed as a functional of the volume 
fraction profiles and SCF potentials for all components in the system. Minimization of this free energy under the 
incompressibility constraint and appropriate boundary conditions yields equilibrium density distributions of various 
components. For most problems of interest, SCF equations need to be solved iteratively and numerically, which 
necessarily involves space discretization, i.e. use of a lattice. Here we employ the method of Scheutjens and Fleer, (lfjj 
which uses the segment diameter a as the size of the cell; throughout this work all the distances are reported in units 
of the cell size a. 

For the present problem, the SCF equations are set up on a cylindrical lattice and solved using two-gradient 
approach, where the density profiles vary both radially and laterally. (l7l - [T9| The lattice layers along the axis of the 
cylinder are numbered according to z = 1, 2, • • • , N z , while circular arrangements of the lattice sites within each z-layer 
are numbered as r = 1, 2, ■ • • , N r . Due to the azimuthal symmetry of the problem, there are L(r) = ir[r 2 — (r — l) 2 ] 
indistinguishable sites at each coordinate (z,r). In order to implement SCF formalism, one needs to define a priori 
step probabilities, which are determined by the fractions of sites in neighboring layers adjacent to a given site on a 
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lattice. In the r-direction, these probabilities follow from geometric considerations: \Vi 



A(r\r-l) = %j-^-, (I) 
3L(r) 

A(r|r + 1) = S ^ 



3L(r)' 

A(r\r) = 1 - A(r|r - 1) - A(r|r + 1), 

where S(r) = 2nr is the surface area between the two cylinders of unit height and radii r and r + 1. The transition 
probabilities along the z-direction are given by:[T3] A(z, z') = 1/3, where z' = z — 1, z, z + 1. Overall, for each lattice 
site (2, r), there are nine transition probabilities given by: 

A(z, r|z', r') = A(z, z')A(r, r') = ±A(r, r'), (2) 

where (z', r') = (z + a,r + ft), with a, /3 = — 1, 0, 1. The transition probabilities thus defined are properly normalized, 
with their sum equal to unity. 

In our model, both lock and key particle are assumed to be uniformly composed of a fixed number of segments 
of type p that fill the corresponding geometrical forms discussed above. Accordingly, the particles are impenetrable 
to any other segment types present in the system, i.e. the volume fraction <p p (z,r)=l for all the lattice sites located 
inside the particles and 4> p (z, r)=0 otherwise. In addition to the particles, there are three other segment types in our 
model, those of free chains (/), grafted chains (g), and solvent (w). The latter fill all the remaining vacancies on the 
lattice, thereby ensuring the incompressibility condition. 

The SCF equation for the volume fraction of grafted chains is written in terms of propagators as follows: fiol Il7l - [l9j 

a (. r A _ r V G g (z,r,s\l)G g (z,r,a\N) 

where the sum runs over all the segments of the grafted chain, the normalization constant C g is determined by the 
grafting density o g , and G g (z,r) is given in terms of the potential u g (z,r): 

G g (z, r) = exp(-f3u g (z, r)), (4) 

where ft = l/k B T. 

The potential u g {z,r) is comprised of two terms: 

u g (z,r) = u hc (z,r)+u g nt (z,r), (5) 

where the hard-core potential u hc (z, r) is independent of the segment type and serves as a Lagrange multiplier enforcing 
the incompressibility condition, while the interaction potential is given by: 

fa g nt (z,r) = ±Y i Xi a <M*,r)>, (6) 

where the sum runs over all the segment types other than g, \ig are the corresponding Flory-Huggins interaction 
parameters, and the step-weighted volume fractions are defined as follows: 

< <t>i( z > r ) >= A ( z ' r ' 2 '' r ')M z ', r ')- ( 7 ) 

z' r' 

As our boundary condition for the grafted chains, we pin their end-segments to the surfaces of lock and key particles 
as depicted in Fig. [TJ [13, GJ, 03 

The SCF equation for the volume fraction of free chains is similar to Eq. ((3]), except that the normalization 
constant is obtained from the value of the bulk volume fraction of free chains (fib (instead of the grafting density), and 
no pinning of the end-segments is performed. Finally, the equation for the solvent is straightforward: 

<j> w (z,r) = exp(-ftu w (z,r)). (8) 



By imposing the incompressibility constraint, we solve SCF equations simultaneously and itcrativcly to obtain the 
equilibrium volume fractions of all the components present in the system. 
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From the resulting volume fractions, one can compute the Helmholtz free energy of the system as follows : [Tol [IH Il9j 

m N, N r ( , , j m \ 

= EEE L ( r )<^> r ) \ ~w + lnG *( z < r ) + E < r ) > f . ( 9 ) 

i — l z—1 r—1 | j>i j 



where the summation index i runs over all the segment types present in the system; Ni = 1, iV, and P for z = w,g, 
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and /, respectively; 4>l is equal to the bulk volume fraction 0& for the free chains, while for the grafted chains 



& = !_Y ( io) 

" E^E^L(r)G g ( Z ,r,N\iy 1 j 

where n is the number of grafted molecules in the system determined by the grafting density. We note that the total 
Helmholtz energy can be split into the contribution due to the grafted chains, f3A g (given by the term i = g in the 
sum over i in Eq. (©), and the contribution due to the mobile species, i.e. free chains and the solvent, j3A w f (given 
by the terms i = w, f in the sum over i in Eq. ([5])). [20] 

The dimensionless lock-key PMF as a function of lock-key separation z is given by the difference of the Helmholtz 
free energy at the particle separation z and at infinite separation. 

/3W(z) = f3A(z) - pA(oc). (11) 

By splitting the free energy into contributions from grafted and mobile species, the PMF can be decomposed into the 
term due to grafted chains and the term due to mobile species: [13, l2p| 

0W(z) = 0W g (z) + 0W mf {z). (12) 

Note that here we focus exclusively on the polymer- and solvent-mediated part of the PMF, i.e. we do not consider 
the bare interaction between lock and key particles, the latter term can always be added separately to the PMF. 



III. RESULTS 



As discussed in the Introduction, prior experimental studies of lock and key colloids have shown that the sharp- 
ness of binding-unbinding transition can be increased by grafting polymeric chains on the surfaces of lock and key 
particles. Furthermore, recent experimental results in the field of polymer nanocomposites indicate that the inter- 
action between sterically stabilized colloids in polymer melts and solutions can be tuned from repulsive to attractive 
by varying the ratio of free and grafted chain lengths. P, [T3T - [l6j Hence, we begin by studying the effect of the ratio 
P/N on the PMF between sterically stabilized lock and key particles. 

The corresponding results are presented in Fig. [3J where we plot dimensionless PMF as a function of the lock-key 
separation for several values of the free polymer volume fraction. In our lattice model, the surface of the bottom of 
the lock cavity corresponds to the surface layer z = 1, and the lock- key separation is measured from the leading edge 
of the key particle, i.e. z = 2 on the graph corresponds to lock and key being in contact. Throughout this study we 
fix the height of the key and the depth of the lock at H = D = 4, the radius of the lock cavity at Ri=6, and the 
length of the grafted chains at N — 10. In the upper panel, we set P/N — 5, and in the lower panel, P/N = 0.5. 
In both cases, the grafting density is fixed at a g = 0.01 and the key particle radius is taken to be Rk = 6, i.e. the 
lock and key sizes match perfectly. One sees that in the case when the length of the free chains is smaller than that 
of the grafted chains, the polymer-induced component is repulsive for all values of the free polymer volume fraction 
considered, which is similar to the behavior of sterically stabilized colloidal particles. [H. I!~3l - fl6| By decomposing the 
total PMF into contributions due to mobile species and grafted polymers (not shown), one finds that the repulsion 
is primarily due to the overlap of the grafted chains, which outweighs the attractive term coming from the osmotic 
pressure of the mobile species in solution (again in analogy to sterically stabilized colloids [l3jV 

On the other hand, when free chains are longer than the grafted ones, the polymer mediated PMF is strongly 
attractive at low polymer volume fraction (06=0.3), moderately attractive at intermediate volume fraction (06=0.7), 
and repulsive at high volume fraction (0b=O.9). We note here that in the case of flat polymer brushes, the phenomenon 
of autophobic dewetting would lead to an attractive PMF in the regime P/N > 1 for all values of free polymer volume 
fractions. [l|, [l3l - [l6j However, for highly curved brushes (grafted on spherical or cylindrical surfaces), melt-brush 
interpenetration and wetting is strongly enhanced. Indeed, recent SCF calculations have shown that wetting-dewetting 
transition is shifted to much larger values of P/N ratio for highly curved sterically stabilized nanoparticles. [2~0] 

This change in the nature of the lock-key PMF from attractive to repulsive with increasing volume fraction 
(when P > N) provides one possible way for tuning interactions in sterically stabilized lock-key colloidal systems. 
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Figure 2 
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FIG. 2. The lock-key PMF between a cylindrical key of radius Rk — 6 and a cylindrical cavity of radius Ri = 6 at three 
different values of the free polymer bulk volume fraction. The length of the grafted chains is N = 10 and their grafting density 
is o g = 0.01. Upper panel: the length of free chains is P = 50, lower panel: P = 5. 



Hence, it would be of interest to decompose the total polymer-induced PMF according to Eq. (|12[) into two separate 
contributions arising due to mobile species and grafted chains. This decomposition is shown in Fig. [31 where the 
contribution due to grafted chains is given in the upper panel, and the contribution due to mobile species in the 
lower panel. Compared to the case when P < N, the situation is reversed in that the term due to grafted chains 
is predominantly attractive, except for a weak repulsive barrier at intermediate lock-key separations at the lowest 
volume fraction considered. The term coming from osmotic pressure due to mobile species is attractive at the low 
volume fraction, weakly repulsive at the intermediate <pb, and strongly repulsive at the highest value of <j>b studied. 
Once again, this behavior is qualitatively similar to sterically stabilized colloids, where it was rationalized by arguing 
that mixing of grafted chains from two particles is more favorable than mixing of short grafted chains and long free 
chains. [l3j 

As a next step, we study the effect of the size matching between the lock and key particles on the PMF. To this 
end, we set the grafting density at a g = 0.05, the length of the free chains at P = 50, and compute the lock-key PMF 
for two sizes of the key particle: Rk — 6 (perfect size matching), and Rk = 3 (key is half the size of the lock cavity). 
Our results are shown in the upper panel of Fig. U for the smaller key, and in the lower panel for the larger key. First 
of all, one immediately notices that the lock-key interaction is substantially weaker in the absence of perfect geometric 
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Figure 3 
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FIG. 3. The decomposition of lock-key PMF shown in the upper panel of Fig. 2 into contributions due to grafted chains (upper 
panel) and mobile species (lower panel). 

matching, in agreement with earlier theoretical and simulation studies. 0, 0, H In addition, while the lock-key PMF 
for the smaller key is weakly attractive at all values of free polymer volume fraction, for the larger key it changes 
from strongly attractive at low volume fraction to repulsive at higher values of <j>b, similarly to what was seen earlier 
in the upper panel of Fig. [5J By comparing the lower panel of Fig. U with the upper panel of Fig. [2J one can see the 
effect of the grafting density on the lock- key PMF (all other model parameters being equal in these two panels). The 
most pronounced difference is the abrupt change in the slope of the PMF at short lock-key separations in the case of 
higher grafting density, which leads to a well-defined minimum in the PMF at the low free polymer volume fraction. 

In all the calculations reported so far, the Flory-Huggins parameters between all the components present in the 
system (free and grafted polymer chains, solvent, and particles) have been set equal to zero, i.e. all the effects observed 
in the lock-key PMFs have been purely entropic. It would be also of interest to consider the role of enthalpic effects in 
tuning lock- key interactions. To this end, we study the effect of specific interactions between free chains and segments 
comprising key and lock particles (by varying the Flory-Huggins interaction parameter between particle segments and 
free chains, Xpf)i an d between free and grafted chains (by varying the Flory-Huggins interaction parameter between 
free and grafted chains, Xgf)- The former results are shown in the lower panel of Fig. [51 and the latter results in the 
upper panel. In these calculations, we have set the free polymer volume fraction at = 0.5, the length of free chains 
at P = 50, the grafting density at a g = 0.05, and the sizes of key particle and lock cavity are perfectly matched. 

All the Flory-Huggins parameters other than those specified in the legend of Fig. O are set equal to zero. Hence, 
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Figure 4 
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FIG. 4. The lock-key PMF for a cylindrical cavity of radius i?;=6 and for two values of the key particle radius: Rk=3 (upper 
panel) and Rk=6 (lower panel). 



the PMF corresponding to Xp/ = i n the lower panel of Fig. [5] can be regarded as that of the reference system in the 
absence of enthalpic effects. This reference PMF displays weak lock-key repulsion at contact and a weak attractive 
minimum at short lock-key separations. In the presence of favorable enthalpic interactions between free chains and 
particles (x P /=-1.5), free chains absorb on particle surfaces, thereby penetrating inside the grafted chains and making 
the lock-key PMF uniformly repulsive. When this interaction is unfavorable (x p /=1.5), free chains are expelled from 
surfaces, which makes the PMF strongly attractive. 



Similar behavior is observed in the upper panel of Fig. [5] where we vary the enthalpic interaction between free 
and grafted chains. When this interaction is favorable (x 9 /=-0.5), the free chains wet the brushes, [l2[ which produces 
a strongly repulsive lock-key PMF. In the opposite case (xgf >0)j the mixing between free and grafted chains is 
unfavorable, which results in a uniformly attractive lock-key PMF, with a deep minimum at short separations. One 
can also note that the range of the PMF in the upper panel is somewhat larger compared to the lower panel, which 
could be due to the changes in the extension of the grafted chains promoted by varying the \gf parameter. 
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Figure 5 
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FIG. 5. The lock-key PMF between a cylindrical key of radius Rk = 6 and a cylindrical cavity of radius Ri = 6. The length of 
the grafted chains is N = 10 and their grafting density is a g — 0.01; the length of free chains is P = 50 and their bulk volume 
fraction is <f>b = 0.5. Upper panel: PMF for three different values of the Flory-Huggins interaction parameter between free and 
grafted chains; lower panel: PMF for three different values of the Flory-Huggins interaction parameter between free chains and 
particles. 



IV. CONCLUSION 

In this work, we have carried out a mean-field theoretical study of lock-key interactions between sterically stabilized 
particles in a polymer solution. Both the key particle and the lock cavity were assumed to have cylindrical shape, and 
their surfaces were uniformly grafted with polymer chains. A two-dimensional version of lattice-based SCF theory 
was employed to compute the PMF between lock and key particles. In analogy with earlier studies of interactions 
between sterically stabilized colloids, it was found that PMF is quite sensitive to several model parameters, making 
lock-key interactions easily tunable. In particular, it was shown that when the free chains are shorter than the grafted 
ones, the PMF is repulsive at all conditions studied, while in the opposite case the lock-key interaction can be tuned 
from attractive to repulsive by increasing the volume fraction of free polymer chains. Likewise, the behavior of the 
PMF, including its range, can be changed by varying the enthalpic interactions between free chains and grafted ones 
or between free chains and particles. Finally, in agreement with earlier experimental and theoretical studies, it was 
observed that the strongest lock-key interaction occurs in the case of perfect size matching between the lock and the 
key; decreasing the size of the key relative to that of the cavity weakens the interaction dramatically. 

It would be of interest to extend the present study to other geometrical shapes of lock and key particles and to 
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go beyond the one-dimensional representation of the PMF by considering different trajectories of the key approaching 
the lock cavity. In addition, it would be of interest to compare the sharpness of binding-unbinding transition (as 
determined by the fraction of occupied lock cavities as a function of free chain concentration) between sterically 
stabilized and charge stabilized systems. This will be the subject of future research. 



V. ACKNOWLEDGMENT 



The author is grateful to Prof. Frans Leermakers for his sharing and help with using of software package sfbox 
which was employed to carry out some of the 2-gradient SCF calculations presented in this work. 



[1] S. K. Kumar and R. Krishnamoorti, Annu. Rev. Chem. Biomol. Eng. 1, 37 (2010). 
[2] S. Sacanna, W. T. M. Irvine, P. M. Chaikin, and D. J. Pine, Nature 464, 575 (2010). 
[3] M. J. Solomon, Nature 464, 496 (2010). 

[4] M. Kinoshita and T. Oguni, Chem. Phys. Lett. 351, 79 (2002). 

[5] N. M. Silvestre, P. Patricio, and M. M. T. da Gama, Phys. Rev. E 69, 061402 (2004). 
[6] P. M. Konig, R. Roth, and S. Dietrich, Europhys. Lett. 84, 5 (2008). 

[7] G. Odriozola, F. Jimenez- Angeles, and M. Lozada-Cassou, J. Chem. Phys. 129, 111101 (2008). 

[8] Z. Jin and J. Z. Wu, J. Phys. Chem. B 115, 1450 (2011). 

[9] S. Asakura and F. Oosawa, J. Chem. Phys. 22, 1255 (1954). 
[10] G. J. Fleer, M. A. C. Stuart, J. M. H. M. Scheutjens, T. Cosgrove, and B. Vincent, Polymers at Interfaces (Chapman and 

Hall, London, 1993). 
[11] A. Striolo and S. A. Egorov, J. Chem. Phys. 126, 014902 (2007). 
[12] A. Milchev, S. A. Egorov, and K. Binder, J. Chem. Phys. 132, 184905 (2010). 
[13] N. Dutta and D. Green, Langmuir 24, 5260 (2008). 
[14] D. L. Green and J. Mewis, Langmuir 22, 9546 (2006). 

[15] A. Bansal, H. C. Yang, C. Z. Li, B. C. Benicewicz, S. K. Kumar, and L. S. Schadler, Nat. Mater. 4, 693 (2005). 
[16] A. Bansal, H. C. Yang, C. Z. Li, B. C. Benicewicz, S. K. Kumar, and L. S. Schadler, J. Poly. Sci. B: Polym. Phys. 44, 
2944 (2006). 

[17] H. P. Huinink, A. de Keizer, F. A. M. Leermakers, and J. Lyklema, Langmuir 13, 6618 (1997). 
[18] B. M. Steels, F. A. M. Leermakers, and C. A. Haynes, J. Chromatography B 743, 31 (2000). 
[19] B. M. Steels, J. Koska, and C. A. Haynes, J. Chromatography B 743, 41 (2000). 

[20] B. van Lent, R. Israels, J. M. H. M. Scheutjens, and G. J. Fleer, J. Coll. Interface Sci. 137, 380 (1990). 
[21] D. M. Trombly and V. Ganesan, J. Chem. Phys. 133, 154904 (2010). 



